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We present a new theoretical approach to the study oiu = 1 quantum Hall bilayer that is based on 
a systematic mapping of the microscopic Hamiltonian to an anisotropic SU(4) spin model on a lattice. 
To study the properties of this model we generalize the Heisenberg model Schwinger boson mean 
field theory (SBMFT) of Arovas and Auerbach to spin models with anisotropy. We calculate the 
temperature dependence of experimentally observable quantities, including the spin magnetization, 
and the differential interlayer capacitance. Our theory represents a substantial improvement over 
the conventional Hartree-Fock picture which neglects quantum and thermal fluctuations, and has 
advantages over long-wavelength effective models that fail to capture important microscopic physics 
at all realistic layer separations. The formalism we develop can be generalized to treat quantum 
Hall bilayers at filling factor u = 2. 
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I. INTRODUCTION 



The physics of high-mobihty two-dimensional electron systems in strong perpendicular magnetic fields continues 
to produce surprises. Electronic states with strong correlations in this regime originate fundamentally from the 
macroscopic degeneracies, Landau levels, in the spectrum of the free-particle kinetic energy Hamiltonian. When spin, 
layer, or other degrees of freedom are important at low energies, correlations at integer filling factors can give rise to 
broken symmetry ground stateall with interesting physical properties. Bilayer quantum Hall systems, whose ground 
states can have spontaneous interlayer phase coherence, provide the simplest and best characterized, example of this 
tendency, and have been studied both theoretically and experimentally, for more than ten yearsB~c3 The superfluid 
effects associated with spontaneous coherence are especially dramatic.E3 Most work on these states has concentrated 
on the case of total lowest Landau level (LLL) filling factor v equal to one. The quantum Hall effect, indicative of a 
gap for charged excitations, is observed in this system at sufficently small interlayer separations d. The charge gap in 
this case must result entirely from interlayer correlations,Ej since individual layer states at filling factor = 1/2 are 
compressible. As the layer separation or the electron density in each layer is decreased there is a phase transition, 
likely of first order in the absejftce-|Of disorder, from a compressible state with no quantum Hall effect to a highly 
correlated incompressible stateoElO with spontaneous coherence. In this paper we present a new approach which can 
be used to study the properties of these interesting ordered states and how they vary with layer separation. 

The layer degree of freedom in the quantum Hall bilayer can be conveniently described using a pseudospin language 
pseudospin up corresponds to an electron in the top layer and pseudospin down to an electron in the bottom layer 
The interlayer phase coherent quantum Hall state corresponds in this language to a pseudospin ferromagnet with an 
X — y easy-plane. Our work is motivated by two limitations that apply to most theoretical work on bilayer quantum 
Hall ferromagnets. The usual assumption in most previous studies of = 1 bilayers has been that the spin degree of 
freedom can be ignored, i.e. that spin dynamics are frozen out by the magnetic field. At finite temperatures, however, 
this assumption is usually not justified: the typical value of the Zeeman gap is only ~ IK, which is often smaller 
than the pseudospin gap when the interlayer phase is fixed by finite tunneling. There is a broad range of nonzero 
temperatures for many samples in which there is a nontrivial interplay between spin and pseudospin fluctuations. A 
second set of limitations that applies to much theoretical work on quantum Hall bilayers follows from the use of a 
gradient expansion in deriving an effective continuum pseudospin model. It is sometimes bothersome in applying this 
approach that divergent coefficients, due to the long-range of the Coulomb interaction, appear in the energy functional 
when one goes beyond the gradient expansion's leading order. A more severe limitation, however, is the fact that this 
approach cannot address the physics that gives rise to the phase transition between ordered and disordered states 
which, as we explain in some detail, is due to competing interactions on microscopic length scales. The importance 
of correlations on microscopic length scales for the phase transition was already evident in early work where it was 
associated with soft collective excitationsB at a finite wavevector. The precise nature of the compressible state close 
to the transition is not completely understood at present and we belive that our approach can shed some light on this 
question. 

In this paper we present a careful study of the spin-pseudospin physics of the v = I bilayer in the quantum Hall 
state which is free from these limitations. Our approach is based on a systematic mapping from the microscopic 
interacting electron problem to a model with SU(4) generalized spins on a lattice and anisotropic interactions. SU(4) 
group appears here since we take into account both spin and pseudospin degrees of freedom and thus single-electron 
states are four-component spinors. The only simplification necessary to derive this mapping is the assumption of the 
absence of charge fluctuations, reasonably well justified in the incompressible state. To take advantage of the existence 
of a charge gap in the quantum Hall regime, we use a complete orthonormal set of magnetic Wannier functions on a 
von Neumann lattice as the LLL orbital basis. Given this basis, and the neglect of charge fiuctuations, we can derive 
a coupled spin-pseudospin Hamiltonian which turns out to be unexpectedly simple. Using the spin-Hamiltonian we 
derive, the full arsenal of analytic and numerical techniques that have been developed to solve spin problems on a lattice 
can be applied to quantuui Hall systems. In this paper, we present a study of this Hamiltonian using the Schwinger 
boson mean field theory,E3 developed originally by Arovas and Auerbach for the SU(2) invariant Heisenberg model, 
appropriately generalized for the anisotropic case. We obtain a variety of new results for the temperature dependence 
of some experimental observables, the spin magnetization and the interlayer differential capacitance in particular. 
We note that our approach is valid not only for the v = 1 case but can be easily generalized to a quantum Hall 
bilayer with any integer filling factor, and with any orbital characteij-Qf nearby Landau levels. For bilayers, the case 
of = 2 is also interesting and has been studied extensively.Eilllj'Eirt3 Since the Schwinger boson mean field theory 
has some limitations, we hope that the present results will motivate the application of more rigorous Monte Carlo or 
exact-diagonalization techniques to our model Hamiltonian. 

The paper is organized as follows. In Sec. || we give a detailed derivation of our spin-pseudospin effective lattice 
Hamiltonian that we obtain by using a functional integral approach. In Sec. Ill we discuss the generalization of the 
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Schwinger boson mean field theory of Arovas and Auerbacfi to the case of spin models with anisotropy. In Sec. IV 
we apply these considerations to our spin-pseudospin model for bilayer quantum Hall systems, discussing the results 
of our calculations and commenting on their relevance to experiments. We conclude in Sec. ^ by briefly discussing 
the application of our approach to other broken symmetry states at integer filling factors in quantum Hall systems, 
including the case of = 2 in bilayers. 

II. GENERALIZED SPIN MODEL 

A single-particle state of an electron in a QH bilayer is described by three quantum numbers: LLL orbital state 
quantum number i, spin a and pseudospin r. In the incompressible state at filling factor v = 1 one can assume at 
low enough temperatures that there are no charge fluctuations and that each LLL orbital state i is always occupied 
by exactly one electron. The model we derive starting from this assumption describes quantum fluctuations of the 
remaining spin and pseudospin degrees of freedom. We will allow electrons in each LLL orbital to be in the most 
general coherent superposition of spin and pseudospin eigenstates. To define our notation, we write the corresponding 
creation operator as 



V'1 = E^^'=4- (1) 



fc=i 



Here we have introduced 4-component spinor index k = 1,...,4 to describe the mixed spin-pseudospin degree of 
freedom, k = 1 labels an up-spin electron in the top layer, k — 2 a, down-spin top-layer electron, fc = 3 a bottom-layer 
up-spin electron, and k = 4 a bottom-layer down-spin electron. c|j, is the creation operator for an electron in the LLL 
orbital state i and with 4-spinor index k. Zik are complex amplitudes that satisfy the normalization constraint 



fc=i 



A single Slater determinant many-body wavefunction can be constructed by putting exactly one electron with an 
arbitrary spinor (|^) (for ly — 1) in every LLL orbital state i 



i*N>-n(E^'^4) io>. 

i \k=i / 



(3) 



Our approach is based on the observation that this set of states is complete when|-eharge fluctuations are neglected. 
We assume that the bilayer is described by a Hamiltonian of the following fornJla 



ki k2i 



+ ^ E Yj 4feicU2^»4fc2C«3fei [{iii2\V+\i3i4) + Tk,ki'rk2k2i^i^2\V-\i3H)] . (4) 

Here we have introduced 4x4 Pauli matrices t° corresponding to the usual Pauli operators acting on pseudospin 
degrees of freedom only. The analogous spin Pauli matrices will be denoted tr". V± = {Vs ± Vd)/2 and Vs and Vd are 
the 2D Coulomb interactions between electrons in the same and different layers respectively. The single particle part 
of the Hamiltonian, h^, consists in general of three terms, corresponding to an interlayer bias potential, an interlayer 
tunneling amplitude, and Zeeman coupling to the perpendicular magnetic field 

^ -^r^ - - (5) 

2 2 2 ^ ' 

Since the occupation of each orbital state is fixed at one for = 1 in our approach, the orbital, or charge degree 
of freedom in is irrelevant. We would therefore like to completely eliminate any reference to the charge degree 
of freedom present in the microscopic Hamiltonian, obtaining an effective Hamiltonian that refers only to spin and 
pseudospin variables. The most convenient and mathematically rigorous way to do this is provided by a functional 
integral approach. Assuming at each discrete imaginary time that the single Slater determinant states are complete. 
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and following standard lines for the derivation of path integral formulations of quantum statistical mechanics problems, 
we obtain the following formal expression for the partition function of a j/ = 1 quantum Hall ferromagnet: 

Z = j D[z,z]e~^^'-'\ (6) 

where the imaginary time action is given by 



Sl'z, z] ^ I dr 



i Y^im, {z,\t^\z,) {z, \t^z,) ~ F+ {z,\z,) {z,\z,) - Fr {z^t^\z,){z, \t^z,)) 



(7) 



Here Hij — {ij\V-\ij) , and Fj^ = (u|^±|j*) Stve the direct and exchange two-particle matrix elements of the Coulomb 
interaction. 

Using this path integral representation as an intermediary, we can now write the Hamiltonian in terms of spin and 
pseudospin operators. The complex fields in this path integral can be thought of as bosonic coherent state variables 
and we can introduce bosonic creation operators a]^^ corresponding to complex variables Zik- The normalization 
condition translates into the single-occupancy constraint for the a bosons 

4 

^alfcOifc^l. (8) 

k=l 

The action is identical to that which arises from a coherent state path integral representation of the partition 
function for a system of bosons described by the following Hamiltonian 



iki k2 



1 



2 /-^ /-^ 

tj ki k2 



[^U«!fci"lfc2«'fci"j'=2^feifei^fc2fe2 - F^^4ky3k2''^k2aJk^ ~ K,^ a.Li "jfe^ ^^^^^ Ajfei ^l^fei ^l,fe,] • (9) 



One can think of the bosons we introduce as Schwinger bosons representing "spins" which are generators of the SU(4) 
group (the usual physical spins are SU(2) generators). The SU(4) generators are written in terms of the Schwinger 
boson creation and annihilation operators as follows 

^■^'^=«k«.fc2- (10) 

The Hamiltonian (9) can now be rewritten in terms of the SU(4) "spin" operators Sf^'^^ 

iki k2 

1 \ ^ \ ^ \ TT cfclfel Cfe2fc2_z _Z TJ+ Qkik2 Qk2ki TP-Qk 

2 [^ y-^i '^'j ^feifei^fe2fe2 ij'^i '^j ij'^t 



2 

ij kik2 



j kiki k2k2 



(11) 



We have thus mapped the original problem of interacting electrons in a quantum Hall bilayer onto a generalized "spin" 
problem. Each LLL orbital i ("lattice site" in the spin model language) is occupied by an SU(4) "spin" and "spins" 
on different sites are coupled via a long-range anisotropic interaction (we will rewrite these spin operators below in 
terms of more physically transparent SU(2) spin and pseudospin operators). 

So far we have not specified the LLL orbital basis we are using. In fact, the accuracy of the model that follows 
from neglecting charge fluctuations does depend on the single-particle representation we employ. It is obvious that 
the usually used orbit-center quantum numbers in Landau or symmetric gauges are not good choices for the complete 
set of orbital labels since they are not localized and the energy penalty for double occupation will vanish in the 
thermodynamic limit. .Qyir. approach is best used in combination with a Wannier-like von Neumann lattice basis 
set for a Landau level,E3Ei in which orbitals are centered on lattice sites and the unit cell area equals the area 
per flux quantum, 2n£'^, to accomodate exactly one electron with a corresponding SU(4) "spin" per site at = 1. 
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We choose a square lattice with lattice constant \J2'kP- . It is not immediately obvious that such a b^is-jexists. A 
strong magnetic field imposes certain restrictions on the localization properties of magnetic orbitalsJ£j~Ej and it is 
well established, for example, that a set of linearly independent ozLijLexponentially localized single-particle orbitals 
in the LLL does not exist. However, it turns out to be possiblet^H to construct a complete orthonormal set of 
Wannier-like eigenfunctions, which, although not exponentially localized, have a well defined Gaussian core and a 
power law fall off at large distances. FoUowingHEj we refer to these orthonormal basis states as magnetic Wannier 



functions. The procedure one uses to construct such a basis set is very much like the one used to construct the usual 
Wannier functions in a crystal. One starts from the set of minimum uncertainty wavepackets for electrons in the 
LLL, centered at the sites of the square lattice descdbed above. The difference from the case of a crystal here is 
that this set is overcomplete, as shown by Perelomov.Ej One then constructs Bloch functions from linear combinations 
of the minimum uncertainty wavepackets and Fourier transforms them to obtain the Wannier functions. Them are 
subtleties in this procedure and we refer the reader to the appendix of our paper and to the original paperg23'E2l for 
further details. For the following discussion it is important to note that the procedure used to construct the basis 
states depends on the size of the system, so that the spin and pseudospin interactions we discuss below have values 
that depend on the overall system size, converging to well defined values in the thermodynamic limit. 

The "spin model" is not very useful as it is because we do not have much intuition about SU(4) "spins". 
However, it turns out to be possible to further pewrite in a much more useful form. To proceed one notes that 
the set of SU(4) generators ( p^ ) is not unique.E3 We define a different set of generators, more appropriate for our 
purpose, as follows. Define total spin and pseudospin operators on each lattice site 



K2 Oifc2 



Ti = ^ a^fc^Tfcjfe^aifcs- (12) 



2 

fclfc2 



It is then possible to show that the original SU(4) generators given by ( p^ ) can be expressed in terms of a new set of 
generators {5°, S'^T''} as follows 



Si' 


= \ + lis-+Tn 


+ sit: 




= l-lis--Tn 


- s!t: 




= l + lis!-Tn 


- S!T^ 


St 


= l-lis!+Tn 


+ StT^ 


sP 


= lst+s+Ti 




sF 






nli 


- StT+ 




c23 

J^ 


- SiT+ 




sr 


= \r+-stTt 




sr 


= \st~StT^^. 





(13) 

Now the SU(4) "spin" Hamiltonian (|ll|) can be rewritten, up to a constant, as the following coupled spin-pseudospin 
Hamiltonian 

H = -Y,[^vn ^ ^tT^ + ^.St] 

% 

I \ ^ (^TJ ^ T?S\rpZrpZ ^ TTi-DrpJ. rp_L ^ Z?-? O O 

ij 

- 2i^f (S, • S,)T^T^ - 2f^^(S, . S,)(T,^ • T|)] . (14) 
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Here F^''° = F+ ± p- . 

This Hamiltonian is one of the main resuhs of our paper. Let us emphasize that this Hamihonian is exact, 
given the assumed absence of charge fluctuations. It is written in terms of the usual physical spin-i operators and 
therefore all the great variety of methods available for spin models can be applied to it. Note also the quite unexpected 



simplicity of (14). All the effects of spin-pseudospin interaction are contained in the two terms with only four spin-^ 
operators. 

The classical ground state of ( p^ ) is one with uniform spin-polarization in the direction of the applied Zeeman field 
and uniform pseudospin that, for Aj = and Ay = 0, is in an arbitrary direction in the x — y plane. It is identical 
to the Hartree-Fock ground state of the original microscopic Hamiltonian. Due to the anisotropic character of ( p^ ) 
there will be corrections to the Hartree-Fock picture even at zero temperature. These corrections can be expected 
to be quite strong, since we are dealing with spin-^ operators, i.e. the system is far away from the semiclassical 
regime. In fact, as we comment later, corrections to Hartree-Fock |theory grow in importance as the layer separation 
is increased. Although these quantum fluctuations can, and have,tJ been addressed in a fully microscopic model at 
the generalized-random-phase-approximation (GRPA) level, the explicit removal of charge degrees of freedom in the 
approximate Hamiltonian we have derived allows a wider variety of more powerful theoretical techniques to be applied. 
In particular, the GRPA corresponds to a linearized spin- wave theory approximation to the above Hamiltonian. Later 
we will go beyond this level approximation, by applying Schwinger boson mean field theory to our model Hamiltonian. 

We now discuss some simple properties that can be used to simplify ( p^ ) before embarking on a more detailed 
study. We note that Ha = , i.e. that on-site direct and exchange Coulomb matrix elements are identical. It 
follows that the on-site contributions to the Hamiltonian add to a constant term that has no effect on spin and 
pseudospin dynamics, a property we use below. We also note that in the absence of the single-particle terms, the 
Hamiltonian has the correct SU{2)spin x U{1) pseudospin symmetry. Correspondingly there are two Goldstone modes 
associated with the spontaneous breaking of these symmetries. Their dispersion relations are readily evaluated in the 
linear spin wave approximation. For Ay = we find that 

Jospin _ A I p+ _ p+ 



i^r""""" = V + + - ^k)' - (Hk - F-y. (15) 

In Fig. |l| we show dispersions ( |l5| ) evaluated for a 20x20 square lattice in the (1, 1) direction using Af = A^ = 0.01 
and d/£ = 1.4. All energies in this paper will be in units of e'^/ei, the characteristic energy scale for all fractional 
quantum Hall systems. Note that in this example, which corresponds to reasonably typical values, the pseudospin gap 
is appreciably larger than the spin gap, even though the bare values of tunneling and Zeeman coupling are the same. 
This difference is due to pseudospin fluctuations present in the ground state, even in a linearized spin-wave theory, 
and originates from the fact that the Hamiltonian is invariant under all spin rotations, but only under pseudospin 
rotations around the z axis, the U(l) symmetry referred to above. The dip in the pseudospin wave dispersion at the 
boundary of the Brillouin zone signals the emerging development of antiferromagnetic instability, which eventually 
destroys the long-range pseudospin ferromagnetic order. It is this feature of the microscopic physics that is missing 
in long-wavelength effective models. 

Let us now briefly consider the simplifled Hamiltonian in which spins are frozen out (S^ • Sj — S^S^ — 1/4), valid 
at T = in the presence of a Zeeman field: 

H^-Y. [AvT^ + A,Tf ] + ^ - F^T^T^ - F^Ti ■ Tf] . (16) 

i ij 

Consider the evolution of (|l^ as one increases the interlayer separation from 0. At d = 0, 2Hij — Ff^ ~ ^F^^ and 
the system in the absence of tunneling and bias is an SU(2) invariant pseudospin fcrromagnet. For < d < d*, 
the F^j term dominates and the bilayer ground state is an easy-plane pseudospin fcrromagnet. At a certain critical 
value of the interlayer separation d* a quantum phase transition occurs to a state with no long-range easy-plane 
ferromagnetic order. To understand more deeply why this occurs we plot in Fig.|^ and ^ the pseudospin couplings 
Jj^^ = 2Hij — Ffj and J^^ = F^ for two different values of the interlayer separation. As can be seen from the 
figures, increasing the interlayer separation primarily affects J^^^ coupling, changing its character from short-range 
ferromagnetic at d/£ = 0.5 to long-range (oc 1/r^) antiferromagnetic at d/£ = 1.4 (we use the word antiferromagnetic 
here just to emphasize the sign change of the interaction, it is not immediately clear that the ground state is actually 
antiferromagnetic, although it very likely is). The change in jf^^ is accompanied by a weakening of the in-plane 
ferromagnetic interaction J^^'^- At present, experimental samples have a d/£ value that is not very far below the 
critical value. We therefore expect that the effective interactions illustrated for the d/£ —1.4 case are representative 
of the current typical experimental situation. This change in character reflects a change in the relative importance of 
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exchange interactions within the layers {Ffj), which are not d dependent and the electrostatic interactions between 
the dipoles {Hij) that are created by z direction pseudospin polarization and whose moments are proportional to d. 
For d/i ~ 1.5, electrostatics dominates, J^^^ is antiferromagnetic, and competition develops between two qualitatively 
different potential ground states. Our analysis makes it clear that it is this developing competition that is responsible 
for the phase transition that occurs in bilayer quantum Hall systems for d/£ ^ 1.5. In our model which neglects 
charge fluctuations, the ground state at large d is likely an antiferromagnet with Ising anisotropy, as it appears 
from the analysis of the classical energetics of Eq.([T^). In experiment, the loss of the quantum Hall effect at this 
transition indicates that it is accompanied by the loss of a charge gap, making charge fluctuations at least somewhat 
important. However, even if we neglect the charge fluctuations, i.e. still assume that the system can be described by 
the pseudospin model ( |l^ in the compressible region, the precise nature of the pseudospin state above the transition 
is not immediately obvious, although as we have pointed out above, it is likely to be antiferromagnetic. The TfTj 
term generally favors an easy axis antiferromagnetic ordering of pseudospins. However, the long-range character of the 
TfTj coupling (since it orginates from Coulomb interactions between dipoles associated with z-direction pseudospin 
polarization, one can expect it to fall off as could introduce frustration and a nontrivial spin liquid state. A 

recent Quantum Monte Carlo studycil of a quantum XXZ model with nearest- and next-nearest-neighbor interactions, 
which is similar to ([l6|), has shown that both Neel and striped phases are possible, depending on the relative strength 
of the interactions. Our model could in principle produce more exotic phases due to the truly long-range character 
ofj-spin-spin interactions in our case. The phase transition between the XY-ferromagnetic and Neel or striped phases 
inE3 was shown to be first order, which gives additional support to the validity of our model, since according to recent 
exact diagonalization studiesJlH the compressible-incompressible transition in = 1 bilayers is very likely of first order 
as well. We believe that the pseudospin Hamiltonian (|l^) captures important parts of the physics of this still poorly 
understood quantum phase transition and deserves further study. 

Although we concentrate in this paper on the case of bilayer quantum Hall systems at = 1 in the lowest Landau 
level, the same formalism can be used to derive spin-pseudospin models for other cases of interest. For example, 
pseudospin stripe stated^ Are expected at = 1 in the case of Landau levels with orbital index N > 2, and Ising 
ferromagnets are expecteda in bilayers when the orbital indices in the two layers differ. In our formalism these 
different phases will appear naturally because of changes in the effective pseudospin-pseudospin interactions. The 
effective interactions we calculate depend slightly on the finite size of the system in which the magnetic Wannier 
states are orthogonalized, but approach a definite value in the thermodynamic limit. Effective interactions for either 
the thermodynamic limit, or for any specified finite size system, are available from the authors on request for any of 
the above cases. 

To investigate finite-temperature properties of our spin-pseudospin model for lowest Landau level bilayers, one 
needs to go beyond linear spin wave theory (|l^) which is not adequate in two space dimensions (2D) at any finite 
temperature due to the absence of long-range order. A simple _h.ut powerful method, the Schwinger boson mean 
field theory (SBMFT), was proposed by Arovas and Auerbach.E£l It does not break spin rotational invariance at 
finite temperatures, yet correctly reproduces linearized spin wave theory in the semiclassical limit. Fo r . th e SU(2) 
Heisenberg model SBMFT qualitatively reproduces the low temperature continuum field theory results.EaO For ID 
integer antiferromagnetic spin chains, SBMFT gives the correct value of the Haldane gap in the ground state. However, 
it fails in the case of the half-integer chains, where Berry phase effects missed by SBMFT result in a gapless excitation 
spectrum at zero temperature. 

In the next section we briefly recount some features of the original theory, developed for the SU(2) invariant 
Heisenberg model, that will be important for our discussion. We then show how to achieve a proper generalization of 
SBMFT for spin models with anisotropy, a procedure that requires some care. We then apply our generalized SBMFT 
to the model Hamiltonian derived in the present section. 



III. SCHWINGER BOSON MEAN FIELD THEORY FOR SPIN MODELS WITH ANISOTROPY 

We first briefly review the original SBMFT proposed by Arovas and AuerbachEl for the SU(2) invariant spin model. 
Our motivation here is to lay the ground work for the generalization to anisotropic interaction models that is required 
for the application to our model Hamiltonian. Consider the 2D Heisenberg model on a square lattice with only 
nearest-neighbor interactions 

F = -ij^S..S,. (17) 

<ij> 

The sum is over neighboring sites with the 1/2 factor correcting for double counting. For J > the ground state 
of this model is ferromagnetic. In 2D long-range order is destroyed at any finite temperature, strong short-range 
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ferromagnetic (or antiferromagnetic for J < 0) correlations are present for temperatures up to ~ | J|, however. A way 
to build a simple theory of this spin liquid state uses the Schwinger boson representation of spin operators 

St = AlO'i2 

a-iflii + a-2ai2 = 25'. (18) 



Correlations in the spin liquid arc then described by bond operators, which are bilinear forms in Schwinger boson 
operators from nearest-neighbor sites. Nonzero expectation values for these bond operators represent short-range 
order. There are of course many possible bond operators one can define, however only the SU(2) invariant ones are 
relevant, since rotational invariance is preserved in the spin liquid state. There are only two SU(2) invariant bond 
operators: 

P}j = Ai^ji + al2«j2 

Aj = 4ia]2 - a\2a]i- (19) 
The operators F"^ and are not independent, but are connected by the following operator identity 

: 4i^, : +4^,, = 45^ (20) 



where colons represent normal ordering. One can write the spin Hamiltonian (|17|) in terms of these bond operators 
as follows 



^ = ~iE H^--4^^.) (21) 



J 

<ij> 



From (|2l|) the physical meaning of the two bond operators is clear. Nonzero expectation value of represents 
short-range ferromagnetic order, of — short-range antiferromagnetic order. Using (20) one can eliminate one of the 
bond operator products from (|2l|), depending on the kind of correlations one expects to have (determined by the sign 
of J in the simple Heisenberg model). Assuming J > one has 



4 ^ 



(: : -252) . (22) 



In the original approach of Arovas and Auerbach this representation is generalized to N Schwinger boson flavors, 
making it possible to set up a systematic expansion of the functional integral representation of the partition function 
in powers of 1/-/V, starting from the saddle point approximation. We are, however, only interested in the saddle 
point, or mean field, solution itself and therefore will use a more simple-minded approach. We assume that the bond 
operator has a nonzero expectation value 

Q^{Fl)^{F,,), (23) 
and perform a Hartree-Fock decoupling of (|2|). The mean field Hamiltonian, neglecting constants, is 



<ji> 



JQ 

<ij> 



^ E ('^Ii^Ji + "I2aj2 + /i.e.) . (24) 



The constraint on Schwinger boson occupation numbers in (|T^) is imposed on average by introducing a chemical 
potential term in the mean field Hamiltonian to yield 

H = (al^aa + ajja^a) - ^ E ("Ii^'ji + '^l2%2 + h.c^ , (25) 

i <ij> 
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or in Fourier space 



= ^ek(akiaki +aj^2"k2), (26) 

k 

where 



ek = A - 2 JQ7k, (27) 

and 7k — ^{cos{kr^a) + cos{kya)) for a square lattice with lattice constant a. It is convenient to redefine A to make 
the notation more physically meaningful, and also to facilitate the solution of SBMFT equations by constraining the 
form of ek 

A^A + 2JQ. (28) 

Then ek becomes 

ek = A + 2J0(l-7k). (29) 
One completes the model by writing down the selfconsistency equations 

^E"^(^k)-2^ 

k 

^ = 25-^E(l-^k)nB(ek), (30) 

k 

where N is the total number of lattice sites. 

The first of equations ( |30| ) is the constraint on the average number of Schwinger bosons per site, while the second is 
the self-consistent expression for the bond operator expectation value. The physical meaning of ek is now clear. The 
second term in ( p9| ) is just the gapless magnon dispersion with a self-consistently renormalized stiffness constant. The 
gap A reflects the finite correlation length in the spin liquid at finite temperatures. As the temperature goes to zero 
and the system size N goes to infinity, A — > 0, Q — s- 25" and the gapless quadratic spin wave dispersion is recovered. 
Note however, that the mean field quasiparticles in SBMFT are not magnons. Due to the fact that the constraint on 
the total number of Schwinger bosons on each site is enforced in (^0|) only on average, unphysical excitations which 
change the local number of Schwinger bosons are allowed in the mean field theory. I— ■ 

It is straightforward to show that at finite temperatures SBMFT reproduces the renormalized classicaEi correlation 
length for 2D quantum ferromagnets and antiferromagnets. For ID integer antiferromagnetic chains at T = 0, the 
Haldane gap is also correcty described. 

Generally, SBMFT has proven to be a very successful and useful tool. One can think of it as an analog of the 
Weiss molecular field theory for quantum spin liquids. It is accurate enough to describe subtle quantum effects, yet 
simple enough to treat complicated spio-jnodels, such as We refer the reader to the original paper of Arovas 

and Auerbach and Auerbach's textbookE3 for further details on the application of SBMFT to the SUli2^ Heisenberg 
model. For a discussion of the application of the SBMFT approach to frustrated antiferromagnets see.E3 

We now generalize the theory outlined above to anisotropic models. It is clear how to do this if one recalls the 
crucial principle of the SBMFT: one needs to write the Hamiltonian in terms of products of bond operators invariant 
under the transformations of the symmetry group of the Hamiltonian. To be concrete, let us consider the following 
model 

E (^'^^'^|-'^^S,^-S|). (31) 

<y> 

This Hamiltonian has the same Z2 x U{1) symmetry as the pseudospin part of our spin-pseudospin model (ph and 
shares the frustration that can occur between antiferromagnetic — interactions and ferromagnetic S — S""*" 
interactions. To illustrate our strategy we assume a square lattice with only nearest neighbor spins interacting. We 
will also discuss only the case J-*- > when the classical ground state is ferromagnetically ordered in the x — y plane, 
which is the case relevant case for bilayer quantum Hall ferromagnets. Other cases can be analyzed in a similar way. 
There are four bond operators that are invariant under symmetry operations of the Hamiltonian: 
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+ al^aj2 




= "I? 




4 




+ al2a]i 


4 







(32) 

The first two are the usual SU(2) ferromagnetic and antiferromagnetic bond operators. The two new bond operators, 
and represent XY and easy axis ferromagnetic correlations. They are connected by the familiar identity 



4j ^ij 



(33) 



As in the SU(2) Heisenberg model case, one can write the Hamiltonian in terms of the products of invariant bond 
operators. There are again several ways to do this, and we need to choose the one that is appropriate to describe the 
correlations in the ground state. In the case where the ground state is an XY-ferromagnet, the desired mean field 
theory follows from the following form for the Hamiltonian 



H 



F^jFij : -XjjXij ] - J- 



<ij> 



(4 



(34) 



We will see later that the resulting Schwinger boson mean field theory correctly reproduces quasiclassical dynamics 
in the S — > 00, T — > limit. 

Regrouping terms and neglecting a constant contribution, the Hamiltonian can be rewritten as 



H = - 



(35) 



To proceed it is convenient to first rotate coordinates by tt/2 around the y-axis. The X'' bond operator then changes 
to 



_ „t „t 



(36) 



F^j does not change of course, because it is an SU(2) invariant. As before, we introduce expectation values of the 



bond operators 

Q = (4) = (F,,) 
P={XI^) = {X,,), 

and perform a Hartree-Fock decoupling. The resulting mean field Hamiltonian is 
H = ^[\-{J^ - J^)Q7k] (aiLiflki + al^ay,2) 

k 

-\{J^ + J"") P^l^ (alio^ki + flkio-ki - al2«-k2 - ak2a- 



(37) 



k2 



(38) 



We introduce the following notation to make the subsequent equations readable 

Aki = A + ( + r-)p + ( - r-)Q{i - 7k) 

Ak2 = A + + ( + r)P + ( - r)Q{l - 7k) 
Skl,2 = ±(J^ + J")P7k 



eki.2 - vAi, 



^kl.2- 



Here we have added magnetic field h in the ai-direction and redefined A as before. 
Eq.(p8|) in this notation becomes 



km"'-kr] 



^km^— k?^ 



k,?n=l,2 



(39) 



(40) 
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Hamiltonian (MG) is diagonalized by a standard Bogoliubov transformation 



cosh 20km = 



cosh6'kmakm + sinheikmalkri 



sinh 29k„ 



Ckn 



and 



"''^km'^km- 



(41) 



(42) 



km 



The parameters in the mean field Hamiltonian are fixed by solving the self-consistency equations 



-y 



^ki 

Ckl 



n.s(eki) + 



1 



Ak2 

ek2 



nB(ek2) + 



1 



-1}=2S 



P 



1 J ^kl 



TV 



n-B(eki) + 



Bk2 

ek2 



Ak2 

ek2 
1 



- 1 



nB(ek2) + 



(43) 



As can be seen from (43), there are corrections to hnearized spin- wave theory even in the T = hmit in the 
anisotropic spin model case. This difference is not unexpected since the XY-ferromagnetic order parameter does not 
commute with the Hamiltonian and has quantum fluctuations in the ground state. If both quasiclassical {S oo) 
and T — > limits are taken, it is easy to see from Eq.(]39|) and (43) that A ^ 0, Q and P 25*, and one recovers the 
magnon dispersion of linear spin wave theory. 

At low temperatures a system described by ( ^l]) must have quasi long range order {i.e. power law correlations) and 
undergo a Kosterlitz-Thouless phase transition at a finite temperature. Unfortunately, our theory fails to reproduce 
these results. Instead our SBMFT always implies a nonzero gap in the Schwinger boson spectrum, and exponential 
correlation functions, at all finite temperatures. This deficiency can be verified by a simple inspection of the first of 
equations (^3|), which represents the constraint on the number of Schwinger bosons, enforced on average. Assuming 
A = 0, the Schwinger bosons dispersion (|3^) is linear for small k. Taking into account only the leading contribution 
to the integral over wavevector, we can take nB(ek) ''^ T/e^. Then it is clear that the integral over k in ( ^ ) diverges 
logarithmically at the origin in 2D. Therefore at any finite temperature the gap in the Schwinger boson spectrum must 
be nonzero, which translates to exponentially decaying spin correlations as the gap introduces a spatial scale into the 
problem. The reason for this failure is the presence of unphysical fluctuations in the local number of Schwinger bosons 
allowed by the soft (enforced only on average) constraint in (^). The existense of quasi long range order depends 
critically on the fact that the physical low energy excitations in ( ^l]) (spin waves) are mostly in the x — y plane. 
Unphysical fluctuations of the Schwinger boson occupation numbers in the mean field theory violate this constraint 
and therefore destroy the quasi long range order. In the SU(2) case, correlations are not so delicate and therefore the 
SBMFT works out qualitatively correctly. 

This deficiency is not a very serious one, however, since in most cases of interest a small pseudospin magnetic field 
is present due to tunneling. We are concerned mainly with calculating quantities like the pseudospin magnetization 
at finite tunneling and the transverse susceptibility, which should not be extremely sensitive to the precise character 
of long-range correlations. We have checked the mean field values of the gap in the XY-model case are indeed much 
smaller at all temperatures than the ones in the isotropic Heisenberg model. As is true for any mean field theory, 
we do not expect our theory to be quantitatively correct, but to capture the general trends in the temperature and 
magnetic field dependence of experimentally observable quantities. Finally, it seems likelythat this problem can be 
corrected by taking into account gaussian fluctuations around the saddle point solution ([43|). Gaussian fluctuations 
of the field A will reintroduce the local constraint and therefore eliminate the unphysical excitations. This issue will 
be examined in the future work. 

We have solved the mean field equations ( p3| ) for finite lattices with up to 100x100 sites and periodic boundary 
conditions which make our grid of wavevectors discrete. Results for the pure XY-model case ( = 0) are summarized 
in Figs. ^ and |^. Fig. ^ shows the in-plane magnetization for different values of the magnetic field in the i-direction. 
Typical temperature dependences of the short-range order parameters Q and P are shown in the inset. At a temper- 
ature of the order of , there is an unphysical transition to the state with no short range correlations. This is a well 
known problem of SBMFT which occurs in the SU(2) case as well. 
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The most interesting quantity for us is the transverse susceptibihty 



as it is the most relevant quantity for experiments in the case of quantum Hall bilayers. It can be written in terms of 
imaginary time spin-spin correlation function as 



^ TV- 



f'^ dT{S!{T)S^^m. (45) 
ij 



The SBMFT result for Eq. (^) can be evaluated using the imaginary time path integral technique. The result is 

-y 

AN ^ 



k 



j4ki ^4k2 ^ -Bki -Bk2 ^ ^ nB{ev2) - ?^s(eki) 



^ki ^k2 , -Bki Sk2 ^ ^^(eki) + '^B(ek2) + 1 



Ckl ek2 Ckl ek2 / ^kl — ek2 

(46) 

Ckl ek2 Ekl ek2 / Ekl + ek2 

The temperature dependence of the transverse susceptibility for the case of XY-model with no in-plane field is shown in 
Fig.^. The monotonic increase of the susceptibility with temperature reflects the softening of the in-plane order due to 
thermal fluctuations. At high temperatures must follow the usual 1/T paramagnetic susceptibility temperature 
dependence. Therefore must have maximum at an intermediate temperatures beyond the range over which 
SBMFT is reliable. Recent quantum Monte Carlo simulation resultscJ support this picture. Our result for thcjsero 
temperature susceptibility = 0.1955 is reasonably close to the value quantum Monte Carlo result 0.2096 fromJiy It 
is interesting to note that in this case SBMFT apparently gives the correct result for the susceptibility, unlike in the 
SU(2) Heisenberg model case_where it is overestimated by a factor of 3/2. The temperature dependence we obtain 
is, however, stronger than ina presumably due to the fact that SBMFT allows unphysical low energy out-of-plane 
spin fluctuations. _ 

In closing this section, we remark that Timm and Jensera have recently suggested a different strategy for generaliz- 
ing SBMFT to anisotropic spin models. These authors employed a more formal 1/7V expansion approach. We believe 
that our approach is better for practical calculations since, unlike the method of Timm and Jensen, it reproduces the 
correct semiclassical dynamics in the limit 5 — ^ oo at the mean field level. 



IV. SCHWINGER BOSON MEAN FIELD THEORY FOR THE v 

SPIN-PSEUDOSPIN MODEL. 



1 BILAYER EFFECTIVE 



In this section we a pply the formalism developed above to the effective spin-pseudospin model of the v = \ quantum 
Hall bilayer given by (|l^. The Hamiltonian (^4|) has an SU(2) invariant spin system coupled to a Z2 x U{\) invariant 
pseudospin system. Therefore our strategy will be to decouple the spin and pseudospin parts of ( p^ in a mean field 
approximation and use SBMFT to study the decoupled spin and pseudospin Hamiltonians. The coupling is restored 
by solving the resulting equations self-consistently. We expect this to be a good approximation due to the strong 
short-range correlations in our model and the long-range character of the couplings in (|lj). 

The decoupled spin and pseudospin Hamiltonians are given by 



2 *J 



2FR{T 



S,; • S9 



(47) 



and 



E 



E 



2H, 



ij^f - 2F^i^, ■ S,)j tit; - Qi^f + 2FPAS, ■ S 



^1 I T," • 



(48) 



To simplify the following equations we define effective couplings for the spin and pseudospin Hamiltonians: 



12 



jf/=2H,,^^Ff^-2Ff^{S,-S,) 



2 

In terms of these couplings (E^ and (E^) are written as 



4 - o^^ + 2^^<7(S.-S,>. (49) 



7J^ = -A.^5f-5]J^S..S„ (50) 



(51) 



We now analyze the effective spin and pseudospin Hamiltonians using the methods developed in section [II. Let 
us consider first. As before, we represent spin operators by two Schwinger bosons, define ferromagnetic bond 
operators and rewrite the Hamiltonian in terms of these operators 



(52) 



The main difference from the simple Heisenberg model we considered before is that here we have long-range spin-spin 
couplings. Correspondingly, bond operators are also defined for all pairs of sites, not just the nearest-neighbor ones. 
We define expectation value of the bond operator for each pair of sites 



and perform a Hartree-Fock decoupling, adding a chemical potential term to account for the constraint 



H 



s 



E 



A 



al^aa + A + 



A, 



(53) 



(54) 



where we have introduced notation J„^ = Jp^Qf^- 

Taking advantage of translational invariance we rewrite (|5J) in Fourier space, simultaneously diagonalizing it and 
reducing the number of variables 



E 



A, 



A Jk aLaki + A + — ^ - Jk aUai^2 



A, 



(55) 



As before, it is convenient to redefine A to make the notation more physical and also to facilitate numerical solution 
of the self-consistent equations. 



A^ A 



(56) 



Then becomes 



(57) 



where 



Ekl — ^ + Jq ~ Jk 

ek2 = A + A^ + - Ji 



(58) 



The two kinds of Schwinger bosons now have different dispersions due to the asymmetry introduced by the magnetic 
field. The Schwinger bosons become the usual magnons in the limit T ^ Q and N ^ oo. On the other hand, the 
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bosons created by a\ Bose condense in the same limit which signals the appearance of long-range order at T = in 
the thermodynamic limit. The copaection between the condensation of Schwinger bosons and long-range magnetic 
order is discussed in more detail inJlS We will always consider finite temperatures and finite system sizes and therefore 
will not be concerned with this. 

The system of self-consistent equations one has to solve now is as follows 



Yl ["s(eki) + ns(ek2)] = 1 

k 

Qk = "s(eki) + ?^s(ek2)- 



(59) 



One therefore has N + 1 coupled equations to solve instead of 2 in the case of the Heisenberg model with nearest 
neighbor interactions. It is useful to define another expectation value which will be used later, 



Qk "B(eki) - "B(ek2)- 



(60) 



These equations have to be solved self-consistently with the corresponding equations for the effective pseudospin 
Hamiltonian that differ from our detailed discussion of the near-neighbor interaction model in a similar manner. The 
effective pseudospin Hamiltonian is 



(61) 



or, in terms of Schwinger bosons and bond operators 



2 ^ 



blibi2 + bl^bii 



jPz) . p.. 



(62) 



Here we use a different notation for Schwinger boson creation and annihilation operators to distinguish them from 
the analogous operators representing real spins. 

It is convenient to first rotate coordinates to diagonalize the tunneling term 



4 ^ 



pt p. . 



SO that the bond operators are 

Introducing expectation values for the bond operators, 



then performing Hartree-Fock decoupling and introducing renormalized couplings. 



J, 



1 



p 



u 2 

we obtain the following mean field Schwinger boson Hamiltonian 

At 



A - 



1 



2 ^ 



[Jij {blibji + bl,b,2 + h.c) + J+ (bib], - blj>l, + h.c 



(63) 



(64) 



(65) 



(66) 



(67) 
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Fourier transforming and redefining A as before 



k 

1 



- -J+ (b^ 

As before, we introduce the notation 



ki - ^k2^-k2 + '^■c- 



Aki = X + J^ + Jo- Jk 
ylk2 = A + At + J+ + Jo" - Jk 
Ski,2 = ±Ji 



(68) 



ekl,2 



^kl,2 ^kl.2- 



Diagonalizing ( |6^ ) by Bogoliubov transformation we finally obtain 

H"" =Y1 (eki/3i^i/3ki + ek2/3i^2/3k2) 



(69) 



(70) 



As usual the parameters in the mean field Hamiltonain are fixed by solving the self-consistency equations 



^ V Uki 



?^i3(ekl) + 2 



Ak2 

ek2 



?lB(ek2) + 2 



-1=1 





Aki 

Ckl 


ns(eki) H 


1' 

" 2 


ek2 


nB(ek2) H 


1" 

' 2 


P^ = 


Ski 

Ckl 


ns(eki) H 


1' 

" 2 


-Bk2 

ek2 


nB(ek2) H 


1" 

' 2 



(71) 



The systems of equations (59) and (|7l| ) have to be supplemented by the equations for the spin-spin and pseudospin- 
pseudospin correlation functions, which couple the spin and pseudospin systems. They can be evaluated straightfor- 
wardly using Wick's theorem. For unbiased bilayers 



P |2 
ij I 



\Q 



x\2 



(72) 



where 



Q 



pp 



{bkib-ki 



+ 6k26-k2' 



^kl 

Ekl 

^ Ski 

Ckl 



nB(eki) - 


1' 




?^s(ek2) ^ 


1" 


^2 




^2 






ek2 





?^B(ekl) + 



Sk2 

ek2 



?T.B(ek2) + 



(73) 



Equations (^9|), (71), ( fz^ ) constitute the full set of equations of the SBMFT for the v — 1 quantum Hall bilayer. 
Fig. 's |6|~[l0| summarize some of the results we have obtained by solving these equations numerically. The calculations 
were performed for lattices with up to 20x20 sites. 

The polarization response of the pseudospin to a field that is transverse to the ordering direction corresponds 
physically to the charge transferred between layers in response to a bias voltage. In the Hartree-Fock approximation 
this quantity has no temperature dependence for temperatures below the charge gap. Our results for the temperature 
dependence of the differential capacitance at d/t = 0.5, 1.0 and 1.4 at several different values of the Zeeman coupling 
strength A^ are summarized in Fig.'s |^-|^. The limit A^ — > 00 corresponds to the usual assumption that the spin 
degree of freedom is frozen out. The large increase of susceptibility compared to the Hartree-Fock value at zero 
temperature in Fig.^ is an artifact of SBMFT, which overestimates the value of the spin-spin correlation function in 
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Eq. ( ^Of ) leading to an overestimate of the transverse susceptibility. We did not try to fix this problem by normalizing 
the value of the correlation function by hand since we are mainly concerned with the qualitative features of the 
temperature dependence of susceptibility rather than exact values for its magnitude. 

One can see that the temperature dependences in Fig.^ are qualitatively different from those in the XY-model 
case m Fig.|. The source of this difference is the dependence of the effective pseudospin interactions in Eq.(^) on 
the spin polarization, as is evident from the strong dependence of the differential capacitance on Zeeman coupling. 
It is clear from thermodynamic considerations that the transverse pseudospin coupling is approximately inversely 
proportional to the difference between the in-plane and transverse effective couplings. From Eq.(49) we see that 
both become less ferromagnetic (or more antiferromagnetic) as correlations in the spin system are suppressed and the 
exchange contributions to the pseudospin effective interactions are reduced. The coupling between spin and pseudospin 
fluctuations therefore tends to make the easy-plane anisotropy strengthen and the transverse susceptibility weaken 
with increasing temperature. This effect competes with the increase in susceptibility with temperature that holds 
for the XY-model with temperature independent coupling constants. Our calculations indicate that the former effect 
dominates as small d/i while the latter effect dominates at larger d/€. This qualitative discussion does not account for 
subtle effects associated with the longer range of the interaction between pseudospin z components, or for role played 
by the relative sign of the XY and Ising pseudospin interaction constants. The transverse pseudospin susceptibility 
can be measured directly by measuring the capacitance of the bilayer system.E3 We emphasize that the temperature 
dependences shown in these figures are all due to correlated quantum and thermal fluctuations and would be absent 
in a Hartree-Fock theory. 

In Fig.^ and |lO|the temperature dependences of the spin and in-plane pseudospin magnetization for different values 
of interlayer separation are plotted. The spin cLeasity in bilayer systems can be measured by measuring the Knight 
shift in an optically pupiped NMR experiments while the pseudospin polarization can be extracted from optical 
absorption experiments. The strong dependence of spin magnetization on the interlayer separation that we find is 
another signature of spin-pseudospin coupling. 



V. CONCLUDING REMARKS 



In this paper we have presented a theory of the finite temperature properties oi v = 1 quantum Hall bilayers 
that accounts for both spin and layer degrees of freedom. The theory is based on an approximate mapping of the 
microscopic interacting electron problem onto a problem of interacting spins and pseudospins on a lattice with a 
prescribed Hamiltonian specified in Eq. (jl^) . Our mapping is exact when fluctuations in the charge density, summed 
over the two layers, are neglected. To investigate the finite temperature properties of this model, we first generalize the 
Schwinger boson mean field theory of Arovas and Auerbach to anisotropic spin Hamiltonians. Using this approach, 
we are able to predict qualitatively the character of the temperature dependence of several experimentally observable 
quantities like the spin magnetization and transverse pseudospin susceptibility, or interlayer differential capacitance. 
The temperature dependences we find reflect a subtle interplay between spin and pseudospin fluctuations. Spin 
fluctuations in bilayer quantum Hall systems can be neglected only at temperatures significantly smaller than the 
unenhanced spin-splitting gap, which is ~ IK in typical circumstances. When spin fluctuations are neglected our 
calculation is closely related to the microscopic generalized random-phase-approximation (RPA) calculation of Joglekar 
et a/.El The linear spin-wave theory approximation to our spin-pseudospin model is equivalent to the microscopic 
generalized RPA; by neglecting charge fluctuations we are able to go beyond this level, approximating the spin-wave 
interactions of bilayer quantum Hall systems using Schwinger boson mean field theory. Our spin-pseudospin model 
is not however coupled in any way to any particular technique and the full array of numerical and analytic methods 
that have been developed for quantum spin models could be employed. 

Our formalism could in principle be generalized to integer filling factors greater than one. To illustrate the com- 
plications involved we discuss briefly the case oi v = 2. Quantum Hall bilayers at = 2 have been studied actively in 
recent ye ars an dpa. number of interesting theoretical results concerning broken symmetry ground states have also been 
obtainedPtMrTllj for this case. In particular, a canted antiferromagnetic phase is predicted in which total spins in the 
two layers have opposite tilts away from the magnetic field and spin-dependent spontaneous coherence is established 
between the opposite layers. Thus spin-pseudospin interplay is important s.t v = 2 even for ground state properties. 
At V — 2 the single Slater determinant wavefunctions that span the Hilbert space of the system have the following 
form 

\m)= n (E<fc4)io), (74) 

i,a=l,2 \fc=l / 
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i.e. there are two occupied orbitals corresponding to each magnetic Wannier state i. In this case, orthonormahty 
conditions for these four orbitals provide four constraints on what will become coherent state labels in the boson path 
integral. Following this line, the 1/ — 2 spin-pseudospin Hamiltonian is 



E 



ij,a0 



- 2i^f (SiQ • ^jp)Tl^TJp - 2F.^{^.,a ■ Sj/5)(T^ • T^^)] . (75) 

Each lattice site now has two spins and two pseudospins corresponding to two occupied orbitals in (Q). It would 
be interesjting to study the correspondence between this Hamiltonian and two similar effective models introduced 
recentlyJlj'CJ Eg. (|75|) is substantially more complicated than ( |l4|) and even the classical ground state, again identical 
to the Hartree-Fock ground state, is known only numerically for general external fields. The presence of implicit 
orthonormahty constraints in Eq.([75|) is a disadvantage of this approach, since it reduces the transparency of the 
Hamiltonian. We believe, nevertheless, that Eq.(|7|) could be be a good starting point for a theory of correlations 
and thermal fluctuations in u = 2 quantum Hall ferromagnets. 
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APPENDIX 



In this section we explain in more detail how our effective interactions are evaluated. A kev_step is the contruction 
of the magnetic Wannier basis, and we start by briefly summarizing the results of Rashba et a/.LJ We want to cons truct 
a complete set of Wannier functions localized at the sites of a square lattice with lattice constant a = v2ttP. We 
start from an overcomplete set of minimum uncertainty wave packets centered on square lattice sites. The minimum 
uncertainty wavepacket localized at the origin is 

cqo is just the zero angular momentum eigenfunction in the symmetric gauge. To construct wavepackets localized at 
other sites of the square lattice we need to translate cqq wavepacket from the origin to each of the lattice sites. Let 
the lattice sites be 



Tmn = ^a,x + nay. (77) 

Define 



yCoo(r) 



(78) 



Here 



^ p-iR(p-fA) 



(79) 



is the magnetic translation operator (this expression is valid in the symmetric gauge only). In the symmetric gauge 



A=|(-2/,x,0). (80) 



Substituting ( |80| ) into (|79| ) one obtains 



17 



^ max nay — V / 

Therefore c„i„ is given by 



exp ( -^rm„ • p ) . (81) 



mn 



c.„(r) = ^=exp 



(r - r^n)^ _j_ 

4^2 + 



z(r X r„ 



(82) 



The functions Cmn are not orthogonal and form an overcomplete set due to the following identity established by 
Perelomov 

oo 

J2 (-l)"+"c™„(r) - 0. (83) 

mn— — oo 

First we will construct a complete set of orthonormal Bloch functions from the set {cm„}. Assume a normalization 
plaquette of area — 2TTi'^N^, where N^f, is the number of magnetic flux quanta. Define a Bloch function at 
quasimomentum k as follows 

^ oo 

*k(r) = ^=== Yl c„„(r)exp(ikr,„„). (84) 



77171 — — oo 



Here J^(k) is a momentum-dependent normalization factor. The allowed values of quasimomentum are determined 
from the boundary conditions. If we assume periodic boundary conditions with respect to magnetic translations 
Tmax,Tnay the allowcd valucs are 

The Bloch functions are assumed to be normahzed to unity over the normahzation plaquette. 

/ dr|*k(r)|' = l. (85) 

J A 

The normalization factor is given by 

oo 

i/(k) = a ^ Cm„(0)cos(kr,„77)- (86) 

mn— — OO 

It turns out that z^(k) goes to zero at the corners of the Brillouin zone krr^ = ±7r/a. At these points the Bloch 
function has to be calculated by a careful limiting procedure (for details see£3). The result is 

oo 

*ko(r) = — 1= i-^r+"^mnc^n{r), (87) 



nin— — oc 



where 7 = —■^J2mn=~ooi~^)"''^^'^mn{0)x'^. Now we can construct magnetic Wannier functions by inverse Fourier 
transform 



(r) = ^= X! *k(i") exp(-zkrm„). 



k 



The functions Wnm are orthonormal and form a complete set by construction. We evaluated our effective spin and 
pseudospin interactions numerically by inserting these functions in direct and exchange Coulomb matrix elements 
{ij\ys,D\'ij) and {ij\Vs^D\ji), where i — {nm). The calculation is best done by Fourier transforming the Coulomb 
interactions Vs.d 

Vs(q = 

27rp2 

^D(q) = ^e-'^ (89) 

and evaluating the plane- wave matrix elements (i|e*'i''|j). The Coulomb interaction matrix elements are then given 
by 
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{ij\Vs,D\ij) = ^ ^(i le'^-"! i) Vs,d{^) {j le-'^^l j), 



(90) 



(91) 
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FIG. 2. Effective pseudospin-pseudospin interaction j/^^ in Eq.jl^) for d/£ = 0.5 (filled circles) and d/£ = 1.4 (open circles). 
Jjj^ changes its character from short-range ferromagnetic at d/£ = 0.5 to long-range antiferromagnetic at d/£ = 1.4, eventually 
making XY-ferromagnetic pseudospin state unstable. 
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FIG. 3. Effective pseudospin-pseudospin interaction in Eq.(p^ j/j'^ for d/l = 0.5 (filled circles) and d/l — 1.4 (open circles). 
Pseudospin XY-ferromagnetic interactions are weakening as the interlayer separation is increased due to the weakening of 
interlayer exchange interactions. 
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FIG. 6. Temperature dependence of the differential capacitance of the bilayer d{T'^) / dAv for d/£ — 0.5 and Az = 0. (filled 
circles), Az — 0.001 (open circles), Az = 0.005 (filled squares) and Az — 0.01 (open squares). The constant Hartree-Fock 
susceptibility is shown by a dashed line. The large increase of the capacitance compared to the Hartree-Fock value is an artifact 
of SBMFT. The main effect determining the temperature dependence of the capacitance at d/£ = 0.5 is the suppression of 
spin-polarization by thermal fluctuations which influences effective pseudospin-pseudospin interactions (see Eq.([49|)). 
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FIG. 7. Temperature dependence of the differential capacitance of the bilayer d{T^) / dAy for d/£ — 1.0 and Az = 0. (filled 
circles), = 0.001 (open circles), A^ = 0.005 (filled squares) and A^ = 0.01 (open squares). The constant Hartree-Fock 
susceptibility is shown by a dashed line. The influence of spin-polarization on the pseudospin system is weaker here compared 
to the case d/£ = 0.5 and the temperature dependence of the differential capacitance is mainly determined by the softening 
of the XY- ferromagnetic pseudospin order by thermal fluctuations, although there still is a substantial dependence on the 
Zeeman-coupling strength. 
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FIG. 8. Temperature dependence of the differential capacitance of the bilayer d{T'^) / dAv for d/£ = 1.4 and Az = 0. (filled 
circles), Az — 0.001 (open circles), Az = 0.005 (filled squares) and Az = 0.01 (open squares). The constant Hartree-Fock 
susceptibility is shown by a dashed line. The spin-polarization influence has become very small both in the temperature and 
Zeeman-coupling dependence. 
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FIG. 9. Temperature dependence of the spin magnetization of the bilayer system for = 0.005, A* = 0.001 and dji = 0.5 
(solid hne), djl = 1.0 (dashed line) and d/£ = 1.4 (long dashed line). The strong dependence of the spin magnetization on the 
interlayer separation is a signature of the spin-pseudospin coupling. 
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